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ABSTRACT 

Recent  improvements  in  rare  earth  magnets  have  made  it 
possible  to  construct  strong,  lightweight,  high  horsepower 
DC  motors.  This  has  occasioned  a  reassessment  of  electrome- 
chanical actuators  as  alternatives  to  comparable  pneuiratic 
and  hydraulic  systems  for  use  in  flight  control  actuators 
for  tactical  missiles.  This  thesis  develops  a  lcw-crder 
mathematical  model  for  the  simulation  and  analysis  of  brush- 
less  CC  motor  performance.  The  model  is  implemented  in  CSMP 
language.  It  is  used  to  predict  such  motor  performance 
curves  as  speed,  current  and  power  versus  torque. 
Electronic  commutation  based  on  Hall  effect  senscr  posi- 
tional feedback  is  sinulated.  Steady  state  motor  behavior 
is  studied  under  both  constant  and  variable  air  gap  flux 
conditions.  The  variable  flux  takes  two  different  forms. 
In  the  first  case,  the  flux  is  varied  as  a  simple  sinusoid. 
In  the  second  case,  the  flux  is  varied  as  the  sum  of  a  sinu- 
soid and  one  of  its  harmonics. 
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I.  INTRODUCTION 

Direct  current  electric  motors  are  gaining  wider  use  ir. 
space  applications  where  long  life,  high  reliability,  high 
torque  and  light  weight  are  critical  factors.  In  recent 
years,  advances  in  magnet  technology  and  materials,  espe- 
cially the  use  of  rare  earth  magnets,  have  made  possible 
significant  improvements  in  the  torque  to  inertia  ratios  of 
permanent  magnet  (PM)  motors.  Samarium  cobalt  is  cue  such 
rare  earth  magnetic  material.  The  great  potential  of  the 
rare  earth  magnets  derive  from  their  inherent  high  flax 
density  and  high  coercivity  properties.  Higher  flux  densi- 
ties mean  greater  developed  motor  torque  while  high  coer- 
civity means  greater  innate  resistance  to  demagnetization 
which  permits  thinner  magnet  sections.  The  two  factors  in 
combination  result  in  increased  mechanical  power  and  energy 
product  with  decreased  physical  size  and  weight. 

Brushless  dc  motors  belong  to  this  class  of  permanent 
magnet  motors  and  thus  enjoy  their  improved  torque  to  size 
characteristics.  In  fact,  brushless  dc  motors  are  signifi- 
cantly smaller  than  either  their  ac  counterparts  or  conven- 
tional brush-type  dc  motors  on  an  equal  horsepower  basis. 
Compared  to  its  conventional  counterpart,  a  brushless  dc 
motor's  structure  is  "inside  out."  In  the  brushless  config- 
uration, the  permanent  magnets  are  attached  to  the  rotcr  and 
the  conducting  coils  are  placed  in  the  stator.  As  its  name 
suggests,  brushes  have  been  eliminated  from  this  type  of 
motor.  Instead,  commutation  of  current  in  the  stationary 
armature  is  accomplished  electronically  by  switching  en  the 
appropriate  windings  via  transistors  specially  designed  to 
provide  the  high  currents  needed  for  motor  power.  In  order 
to  sequence  the  current  to  the  coils,  the  commutation 
circuitry  is   provided  rotor  position  feedback   from  sensing 


devices,  most  commonly  from  Hall  effect  sensors.  These  are 
simple  semiconductor  devices  that  develop  a  pclari2ed 
voltage  depending  upon  the  magnetic  field  passing  through 
them.  Small,  highly  sensitive  and  reliable.  Hall  effect 
senscrs  require  little  power  to  operate  and  thus  are  widely 
used. 

Brushless  dc  motors  possess  all  the  advantages  tradi- 
tionally associated  *ith  conventional  dc  motors  such  as 
linear  torque-speed  characteristics,  excellent  response 
times  and  high  efficiency.  In  addition,  there  are  a  number 
of  important  advantages  that  brushless  dc  motors  possess 
over  standard  brush-type  dc  motors.  An  obvious  advantage  is 
that  electronic  commutation  means  the  elimination  of  commu- 
tators, the  wear  and  tear  of  brushes,  and  the  build  up  of 
carbon  dust  which  heretofore  significantly  limited  the  life 
of  the  standard  dc  motor.  Since  electronic  currenb 
switching  essentially  involves  no  component  wear,  there 
results  improved  reliability  and  enhanced  lifespan.  As  well 
as  being  subject  to  wear  and  tear,  the  action  of  brushes 
sliding  over  commutators  also  generates  radio  frequency 
interference  (RFI)  which  is  a  severe  liability  in  many 
applications.  In  fact,  it  is  often  in  hostile  and  explosive 
environments  that  the  advantages  of  dc  motors  are  most 
needed.  Another  benefit  is  in  heat  transfer  characteris- 
tics. By  placing  windings  in  the  stator  the  thermal  path  to 
the  environment  is  made  shorter  and  more  direct.  With  the 
removal  of  heat  thus  enhanced,  mechanical  and  electrical 
malfunctions  due  to  heat  are  reduced.  A  further  important 
advantage  of  the  brushless  dc  motor  is  in  the  area  of  motor 
speed  variability.  Whereas  ac  motors  usually  operate  at  one 
speed  fixed  by  the  powerline  frequency,  any  dc  motor's  speed 
can  be  changed  to  meet  varying  power  requirements  by  simply 
adjusting  the  dc  input  voltage.  The  electronic  commutation 
scheme  of  the  brushless  dc  motor  lends  itself  nicely  to  just 


10 


this  type  of  variable  speed  control.  It  is  a  relatively 
simple  matter  to  attach  to  the  commutation  circuit  beard  a 
small,  inexpensive  microprocessor  whose  logic  can  be 
programmed  to  run  the  motor  at  any  of  several  speeds  or  in 
either  direction.  A  single  motor  can  thus  be  programmed,  in 
place  if  desired,  so  that  its  speed,  direction  and  tcrgue 
are  matched  to  a  particular  application. 

There  are  at  least  as  many  applications  for  brushless 
dc  motors  as  there  are  for  brush-type  motor  systems  as  well 
as  a  hest  of  new  applications.  Of  particular  relevance  to 
this  thesis  is  the  use  of  brushless  dc  motors  in  flight 
control  actuator  sytems  of  Navy  tactical  missiles,  in 
particular,  the  cruise  missile.  There  are  three  general 
types  of  actuators  that  are  in  use  in  missile  systems:  pneu- 
matic, hydraulic  and  electromechanical.  The  maneuvering 
requirements  of  the  mere  advanced  tactical  missiles  call  for 
high  control  torques  which  heretofore  could  only  be  produced 
by  high  pressure  pneumatic  or  hydraulic  actuators.  3cth  of 
these  actuator  systems  are  more  susceptible  to  mechanical 
malfunction  than  their  electrical  counterparts.  Kith  the 
improved  torque  to  inertia  ratios  resulting  from  better  rare 
earth  magnet  materials  and  the  higher  mecnanicai  reliability 
due  to  removal  of  brushes,  as  well  as  other  advantages 
previously  enumerated,  the  brushless  dc  motor  is  in  a  posi- 
tion to  compete  with  most  pneumatic  and  hydraulic  systems. 
This  thesis  is  but  one  part  of  a  detailed  study  of  the  state 
of  the  art  in  electromechanical  actuators  as  applied  tc  the 
cruise  missile  system  [Bef.  1 ]. 


1  1 


II.  NATURE  OF  THE  PROBLEM 

A.   PROGRAMMING  LANGUAGE 

The  objective  of  this  thesis  is  to  create  a  simplified 
low-order  mathematical  model  that  will  accurately  simulate 
the-  response  of  a  brushless  dc  motor.  The  IBM  Continuous 
System  Modeling  Program  (CSMP)  was  chosen  for  use  since  it 
is  a  user-oriented  computer  language  specifically  created 
for  modeling  dynamic  physical  systems  and  is  one  of  the  most 
widely  used  in  this  area  £Ref.  2].  It  has  the  flexibility 
of  determining  the  response  of  a  system  that  is  modeled 
either  in  a  block  diagram  format  or  as  a  set  of  ordinary 
differential  equations.  While  CSMP  is  an  offshoot  of  IBM's 
Digital  Simulation  language,  it  embodies  Fortran  as  its 
source  language  and  retains  much  of  Fortran's  capability  and 
flexibility.  CSMP  is  an  applications-oriented  program 
intended  for  use  with  the  IBM  mainframe  systems.  The 
utility  and  ease  of  using  this  program  language  stems  from 
simplified  program  control  statements  that  almost  exactly 
describe  the  mathematical  equations  or  physical  variables  of 
the  system  and  from  the  flexibility  of  its  program  structure 
which  contains  preprogrammed  function  blocks  that  obviate 
the  need  for  complicated  subroutine  programming.  In 
essence,  CSMP  was  chosen  because  it  permits  the  user  to 
concentrate  on  the  details  of  the  physical  system  rather 
than  on  the  time-consuming  complexities  of  programming. 

E.   SISTEM  BLOCK  DIAGBAM 

The  block  diagram  format  was  chosen  as  the  basis  for 
the  CSMP  program  since  it  is  a  convenient  way  of  analyzing 
the  input-output  relationships  of  a   physical  system  and  the 
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flow  of  signals  through  it  without  having  to  refer  to  the 
system  equations.  The  basic  idea  of  a  block  diagram  stems 
from  the  application  of  Laplace  transforms  to  a  system's 
integro-diff erential  cr  differential  equations.  In  effect, 
a  differential  equation  is  rendered  into  a  simpler  algebraic 
expression  which,  in  turn,  becomes  a  fundamental  equation  of 
a  system's  functional  blocks.  Each  block  indicates  the 
"relationship  between  its  input  driving  signal  and  output 
response;  ie.,  the  transfer  function  of  the  input  and 
output.  The  transfer  function  is  defined  as  the  ratio  of 
the  Laplace  transform  of  the  output  to  the  transform  of  the 
input  with  all  initial  conditions  assumed  to  be  zerc.  A 
block's  transfer  function  does  not  include  any  information 
about  the  internal  structure  of  that  part  of  a  system, 
merely  the  transformation  of  a  signal  between  two  points  in 
a  system.  Furthermore,  it  should  be  borne  in  mind  that  all 
physical  systems  have  certain  transfer  characteristics  that 
are  actually  nonlinear  to  some  extent.  A  ic  motor  system 
can  be  fairly  accurately  described  using  the  transfer  func- 
tion approach  to  modeling  and,  when  nonlinear  transfer  char- 
acteristics exist,  CSMP  modeling  can  usually  be  obtained 
from  preprogrammed  nonlinear  function  routines. 

C.   DEVELOPMENT  OF  HCTOR  SYSTEM  EQUATIONS 

In  order  to  construct  the  block  diagram  of  a  system, 
one  approach  is  to  identify  the  various  mathematical  rela- 
tionships between  the  components  of  the  system  and  to  write 
balance  eguations  that  will  identify  each  individual  block 
comprising  the  system.  The  system  equations  for  a  brushless 
dc  motor  are  presented  in  chapter  three  where  ail  the 
details  of  the  logic  control  and  switching  of  transistors 
for  commutation  of  the  motor  are  developed.  What  follows  is 
a  simplified  input-output  characteristics   model  of  the  type 
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that  has  been  a  standard  for  brush-type  dc  motors  used  in 
control  system  studies  [Ref.  3].  The  development  of  this 
input-output  model  is  the  first  learning  step  in  the 
modeling  process  for  the  brushless  dc  motor  and  has  the 
additional  usefulness  of  being  a  tool  for  a  concurrent 
mechanical  engineering  study  cf  load  torque  requirements  for 
the  motor  [ Ref .  1 ]. 

To  write  the  electrical  balance  equation  describing  a 
basic  dc  motor,  the  classical  loop  method  based  on 
Kirchoff's  voltage  law  is  applied  to  the  equivalent  motor 
circuit  (figure  2.1)  . 


Figure  2. 1    Equivalent  DC  Motor  Circuit 
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The  basic  dc  motor  has  the  coil  windings  on  the  armature  and 
the  magnetic  field  source  on  the  stator.  In  the  case  of  the 
brushless  dc  motor,  the  rotor  consists  of  permanent  magnets 
that  produce  the  magnetic  field  and  the  stator  contains  the 
coil  windings.   The  Kirchoff  Law  for  the  stator  circuit  is 

es  (t)  =  Ldi/dt  +  Ri(t)  +  eb(t)  (eqn  2.1) 

where  R  and  L  are  the  stator  resistance  and  inductance 
respectively,  and  eb  (t)  is  the  back  emf.  When  a  conductor 
moves  in  a  magnetic  field,  or  there  is  any  relative  action 
between  the  magnetic  field  source  and  the  conductor,  a 
voltage  is  generated  across  the  terminals  of  the  conductor. 
In  the  case  of  the  dc  motor,  the  voltage  is  proportional  to 
the  shaft  velocity  and  tends  to  oppose  the  current  flow. 
The  relationship  between  the  tack  emf  and  the  shaft  velocity 
is 

es(t)  =  (Km*(J))  *wm  (t)  =  Kb*wm  (t)  (eq-n  2.2) 

es(t)  =  Kb*d0m  (t)  /dt  (eqn  2.3) 

Substituting  and  rearranging  equations  2.1  and  2.2  into 
their  differential  forms  gives  the  following  two  equations. 

di/dt=  (1/L)*es(t)  -  (E/L)*i  (t)-(1/L)  *Kb*wm  (t)      (eqn  2.4) 

dem(t)/dt  =  wm(t)  (eqn  2.5) 

The  basic  balance  equation  that  describes  the  motor  as 
a  mechanical  system  derives  from  application  of  Newton's 
laws  of  motion  to  the  system  components.  The  dynamic  equa- 
tion   for   a    motor  coupled    to    its    load    is 

tm(t)     =    J*dwm(t)/dt    +    B*wra  (t)     +    tl(t)  (eqn    2.6) 

15 


where  tl  (t)  is  the  load  torque.  Rearranging  the  equation  it 
becomes 

dwm  (t)/dt=(1/J)*tm  (t)-(1/J)  *tl  (t)  -  (B/J)  *wffi  (t)        (eqn    2.7) 

E  is  the  total  viscous  friction  coefficient  of  the  motor  and 
is  comprised  of  the  sum  of  the  viscous  friction  due  to  the 
motor,  Bm,  and  the  load,  Bl,  as  seen  by  the  motor  shaft. 
Likewise,  J  is  the  total  system  inertia  comprised  of  the  sum 
of  the  motor  inertia,  Jm,  and  of  the  load  inertia,  Jl, 
reflected  through  the  coupling  device  to  the  motor  shaft. 
Eoth  Bl  and  Jl  are  functions  of  the  motor's  speed  reducing 
mechanism  where  N  is  the  speed  reduction  ratio  (K>1).  The 
following  equations  summarize  the  foregoing. 

3  =  Bm  +  Bl  (eqn  2.8) 

J  =  Jm  +  Jl  {ein  2.9) 

Bl  =  Blp/N  (eqn  2.  10) 

Jl  =  Jlp/N  (eqn  2.  11) 

Jlp  and  Elp  are  the  inertia  and  viscous  friction  seen  at  the 
load  side  of  the  system.  Since  a  dc  motor  is  essentially  a 
torque  transducer  that  converts  electrical  energy  into 
mechanical  energy,  the  following  equation  is  necessary  to 
show  the  relationship  amoung  the  developed  torque,  tm(t), 
the  air  gap  flux,  $,   and  the  stator  current,  i (t) . 

tm(t)  =  (Km*<J))  *i(t)  (eqn  2.12) 

Because  the  magnetic  field  in  the  motor  is  assumed  uniform 
and  constant,  this  can  be  simplified  to 
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tm  (t)  =  Kt*i  (t)  {egn  2.  13) 

vhere  Kt  is  the  torgue  constant  of  the  motor. 

Equations  2.1  through  2.13  represent  the  cause  and 
effect  equations  of  the  system.  The  application  of  es  (t) 
produces  a  current  flow  which  causes  the  torgue  and  the  back 
emf  tc  b€  generated.  The  torgue  produced  then  causes  the 
angular  displacement  6m  (t) .  Given  the  differential  equa- 
tions of  the  system,  the  Laplace  transform  can  be  applied 
and  the  block  diagram  of  the  system  generated.  As  an 
example,  the  block  for  eguation  2.1  is  generated  as  follows 

es  (t)  -  eb(t)  =  L*di/dt  +  R*i{t)  (egn  2.14) 

^{es(t)  -  eb(t)}  =XlL*di/dt  +  R*i  (t)  }  (egn  2.15) 

En  (s)  =  Es{s)  -  Eb  (s)  =  (Ls  +  R)  *I  (s)  (egn  2.16) 

I(s)/En(s)  =  Is/(Es  -  Eb)  =  1/(Ls  +  R)  (egn  2.17) 

Performing  the  same  operation  to  the  remaining  differential 
eguations,  the  block  diagram  in  Figure  2.2  results. 

It  is  important  to  note  that  while  a  dc  motor  is  basi- 
cally an  open-loop  system,  the  back  emf  of  this  type  of  dc 
motor  acts  as  a  natural  feedback  loop  that  tends  to  improve 
the  stability  of  the  motor. 

Cnce  the  basic  block  diagram  is  developed,  writing  a 
CSMP  program  is  simple  and  fairly  straightforward.  Care 
must  be  taken,  however,  to  put  the  block  transfer  functions 
into  the  proper  format  for  CSMP's  functional  blocks.  A  copy 
of  the  CSMP  program  for  a  basic  dc  motor  is  provided  in 
Appendix  A. 
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Figure  2. 2   Block  Diagram  of  DC  Motor 


III.  CORYE  PREDICTION 

A.   OVERVIEW 

Cue  of  the  purposes  of  this  model  is  to  produce  a  set 
of  performance  curves  that  accurately  simulate  a  given 
motor's  behavior  when  operated  under  varying  load  condi- 
tions. As  such  the  iiodel  can  be  a  useful  tool  for  studying 
the  changing  performance  due  to  the  configuration  changes  of 
a  motor  under  development  or  for  observing  the  behavior  of  a 
motor  under  varying  conditions  of  application  or  environ- 
ment. Thus,  having  developed  an  input-output  computer  model 
for  the  EC  motor,  the  next  step  was  to  assign  values  tc  the 
model  parameters  and  run  the  program.  Program  inputs, 
constants  and  parameters  were  assigned  values  that  are 
typical  of  brushless  DC  motors  commercially  available. 
Table  I  provides  the  parameter  values  for  a  given  motor  that 
will  te  used  as  a  generic  motor  for  the  following  analysis 
and  discussion. 


TABLE  I 
Typical  Commercial  Motor  Parameters 


Stator  Resistance,  R 
Stator  Inductance,  L 
Torque  constant,  Kt 
Back  EHF  constant,  Kb 
Rotor  Inertia,  Jm 
Viscous  Friction 

Coefficient,  Em 


2.74 

ohms 

0.0016 

henries 

15.9 

oz-in/amp 

0.  112 

volt/rad/s 

0.001 

oz-in/s2 

unknown 

oz-in/rad/ 

19 


B.       SPEED    VS    TOEQUE    CURVE 

Specifically,  three  curves  were  produced:  motor  speed 
vs.  generated  torque,  current  vs.  torque  and  output  power 
vs.  torque.  Since  permanent  magnet  DC  motors  have  linear 
characteristics,  a  straight- line  curve  resulted  under  all 
load  conditions  when  the  parameters  of  Table  I  were  inserted 
into  the  model  and  a  voltage  of  30  VDC  was  applied.  In 
order  to  fully  understand  the  individual  contributions  of 
the  various  motor  parameters  on  the  speed-torque  curve,  it 
was  necessary  to  determine  the  mathematical  relationships 
between  the  parameters.  Referring  to  equations  2.1  and  2.5, 
the  electrical  and  motor  equations  when  considered  under 
steady   state  conditions    become 

es  (t)  =   E*i(t)    +    Kb*wm(t)  (eqn    3.1) 

tm (t)  =   Kt*i  (t)  (eqn   3.2) 

Solving  for  i{t)  and  combining  the  equations  results  in 

es  (t)  =  (R/Kt)*tra(t)  +  Kb*wm(t)  (eqn  3.3) 

When  there  is  no  torque  applied  to  the  motor,  the  current  is 
small  enough  to  neglect  and  the  no-load  velocity  becomes 

wnl  =  es/Kb  (eqn  3.4) 

Given  that  the  applied  voltage  is  a  constant,  it  is  apparent 
that  knowledge  of  the  value  of  the  back  emf  constant,  Kb, 
will  provide  the  various  no-load  speeds  for  any  value  of 
applied  voltage.  On  the  other  hand,  in  the  absence  of  any 
known  motor  specifications,  the  no-load  speed  of  a  given 
motor  can  be  empirically  measured  and  then  the  back  emf 
constant,    Kb,    can   be    mathematically    derived. 
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Looking  next  at  the  slope  of  the  speed-torque  curve,  it 
can  be  seen  that  whereas  one  endpoint  of  the  curve  termi- 
nates at  no-load  speed,  the  other  ends  at  stall  which  is 
that  value  of  torque  which  is  large  enough  to  literally  stop 
the  motor.  From  equations  3.1  and  3.2,  with  wm  set  to  zero 
the    generated    torque,    tm(t),    at    stall   is   ts   given    by 


ts    =    (Kt/R)*es 


(eqn    3.5) 


Combining   equations    3.4   and    3.5,       the   equation   of    a   straight 
line    for   the   speed-tcrque   curve    results    in 


wm  (t)     =   wnl   -     (Rm*toa(t)) 


(eqn    3.6) 


where 


Rm    =   E/(Kb*Kt) 


(eqn    3.7) 


and  Em  is  the  slope  factor  for  the  curve,  also  called  the 
speed  regulation  constant.  Obviously,  changing  any  or  all 
of  the  parameters  will  alter  the  slope  factor,  Rm.  However, 
the  back  eif  constant,  Kb,  and  the  torque  constant,  Kt,  are 
strictly  related  by  the  common  factor,  air  gap  flux,  (J).  In 
fact,  when  both  constants  are  put  in  the  MKS  system  of 
units,  Kb  in  volt/rad/s  and  Kt  in  Nm/arap,  then  they  are 
equal  as  equation  3.8  shows. 


Kb  (volt/rad/s)  =  Kt  (Nm/amp) 


(eqn  3.8) 


Thus,  knowledge  of  either  constant  provides  knowledge  of  the 
other.  Referring  to  equation  3.7,  it  can  be  seen  that  the 
final  unknowns  are  the  resistance,  R,  and  the  slope  factor, 
Rm.  If  the  resistance,  R,  is  known  from  manufacturer  speci- 
fication sheets,   for  example,   then  the  slope,   Rm,   can  be 


21 


easily  solved.  On  the  other  hand,  if  given  an  empirically 
measured  speed-torque  curve,  it  is  now  obvious  that  the 
curve  not  only  provides  Kb  which  in  turn  provides  Kt,  but  it 
also  can  provide  a  derived  value  for  motor  resistance,  R. 
Once  acquired,  the  three  basic  parameters,  Kb,  Kt  and  R, 
were  plugged  into  the  CSMP  model  and  the  speed-torque  curve 
f or • a   30    VDC   input    was   generated    (see   Figure    3.1). 

This  same  procedure  can  be  used  to  produce  a  family  of 
speed-torque  curves  if  it  is  desired  to  study  a  motor's 
performance  under  different  input  voltages.  Figure  3.2 
shows  the  family  of  curves  produced  when  the  values  given  in 
Table    I   were   entered    into    the    model. 

C.  CURRENT  VS  TORQUE  CURVE 

With  the  speed-tcrque  curve  understood,  the  next  task 
was  to  repeat  the  same  general  procedure  for  the  motor 
current  vs.  torque.  As  a  starting  point,  the  equation  for 
the   steady-state  current   is   given   by 

iss  (t)    =    (es(t)    -    Kt*wm(t))/R  (egn    3.9) 

Since  Kt  and  R  were  fixed  by  the  speed- torque  curve,  gener- 
ating a  current-tor gue  curve  was  straightforward  procedure 
for    any    applied    voltage    (see    Figure    3.3). 

D.  OUTPUT    POWER    VS    TCRQDE   CURVE 

The  final  performance  curve  to  be  studied  was  the 
output  pcwer  vs.  torgue  curve.  The  basic  equation  fcr  power 
is  given  by 

P(t)  =  tm(t)*wm(t)  (eqn  3.10) 

P(t)  =  Kt*i(t)  *wm(t)  (eqn  3.11) 
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Figure  3.1    Motor  Speed  vs  Torque  Curve 

where  pcwer  is  measured  in  watts.  Since  this  equation 
involves  do  new  constants,  producing  the  power  curve  was 
only  a  matter  of  adding  a  line  to  the  existing  model  which 
also  included  a  conversion  factor  for  changing  the  torque 
units  frcm  oz-in/amp  to  Nm/amp  in   order  for  the  power  to  be 


23 


SPEED-TORQUE  CURVES 


TORQUE  (OZ-IN) 


192 


Figure  3.2    Family  of  Speed-Torque  Curves 

dimensionally  correct  in  watts.  Again  running  the  program 
with  the  parameters  given  in  Table  I,  the  model  produced  the 
power  curve  given  in  Figure  3.4 
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CURRENT  VS  TORQUE  CURVE 
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Figure  3.3    Motor  Current  vs  Torque  Curve 

As  expected  both  the  speed-torque  and  cue  en  t-  torque 
curves  were  linear.  The  power-torque  curve  is  obviously  not 
linear,  especially  at  high  leads  where  the  power  begins  to 
roll  off  fairly  sharply.  From  the  modelling  point  of  view 
this  is  due  solely  tc  the  fact  that  the  motor  slews  dewn 
faster  than   the  load  torque   increases  (see   equation  3.9) . 
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Figure  3.4    Output  Power  vs  Torque  Curve 

Another  factor  that  further  accentuates  this  power  rollcff 
in  an  actual  motor  is  the  effect  due  to  armature  reaction. 
Whenever  a  current  flows  through  the  motor  armature  in  a 
permanent  magnet  dc  motor,  the  armature  becomes  an  electro- 
magnet which  produces  a  flux  that  tends  to  oppose  the  flux 
produced  by  the  permanent  magnets.     This  has  the  effect  of 
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partially  demagnetizing  the  permanent  magnets.  This 
demagnetization  is  a  reversitle  effect,  meaning  that  as 
current  returns  to  zero  the  permanent  magnets  return  to 
their  full  strength.  However,  as  the  current  increases  and 
armature  reaction  sets  in,  it  has  the  effect  of  decreasing 
the  torgue  constant,  Kt  (as  well  as  Kb).  As  equation  3.9 
indicates,  this  also  contributes  to  the  output  power  roiicff 
at  high  loads.  Though  current  permanent  magnet  DC  motors 
made  of  rare  earth  magnets  have  high  coercivity  that  resists 
armature  reaction,  at  high  levels  of  rated  current  a  small 
amount  of  armature  reaction  still  occurs.  For  the  purposes 
of  this  thesis,  the  irodel  is  considered  to  accurately  simu- 
late motor  behavior  for  loads  up  to  the  peak  power  load  tor 
the  given  applied  voltage.  From  that  load  upward,  it  is 
assumed  that  armature  reaction  is  likely  to  occur  and  result 
in  nonlinear  behavior  which  is  not  included  in  this  model. 

E.   MCTOB  HEVERSA1 

The  final  step  in  the  modeling  procedure  was  to  reverse 
the  motor's  direction  and  ensure  that  mirror  images  of  the 
three  curves  resulted.  To  do  this,  it  was  first  necessary 
to  replace  the  positive  input  voltage  with  a  negative  count- 
erpart. Next,  in  order  to  maintain  model  consistercy  where 
the  lead  torque,  tl(t),  is  treated  as  a  torque  that  opposes 
the  motion  of  the  mctor,  it  was  necessary  to  codify  the 
model  so  that  the  load  torque  was  a  positive  value.  This 
consisted  of  creating  a  CSNP  procedure  block  that  identified 
the  input  voltage  as  either  positive  or  negative  and  then 
treated  the  load  torgue  accordingly.  Once  program  tugs  were 
removed,  the  model  accurately  simulated  motor  reversal  with 
proper  speed,  current  and  power  values  for  the  appropriate 
load  conditions. 
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IV.  ELECTRONIC  COMMUTATION 

A.   SWITCHING  ACTION 

The  next  major  step  in  this  thesis  was  to  convert  the 
basic  model  of  a  standard  brush-type  dc  motor  into  a  three- 
phase,  four-pole  brushless  dc  motor  system.  Figure  4.1 
shows  this  system  with  the  three-phase  stator  windings 
configured  in  a  standard  'star'  connection  with  each  winding 
oriented  120  electrical  degrees  from  the  others. 


Figure  4.  1   Controller  Configuration 


The  six  transistors  are  connected  to  the  ends  of  each  stator 
leg  and  through  logic-controlled  switching  action  provide 
three-phase  full-wave  motor  control.  Figure  4.2  indicates 
the  switching  logic  sequence  for  counterclockwise  rotation. 

To  produce  clockwise  rotation,  the  sequence  is  reversed. 
In  either  case,  a  pair  of  transistors  will  be  switched  at 
the  start  of  each  30  degree  (mechanical)  interval  which 
causes  current  to  simultaneously  build   up  in  one  leg,   flow 
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Figure  4.2   Switching  Logic  of  a  3-Phase  Brushless  DC  Motor 

steadily  in  a  second,  and  decay  to  zero  in  a  third-  Given 
proper  sequencing,  the  net  developed  torque  ideally  reaches 
a  steady  state  value  that  causes  motor  rotation  at  a 
constant  speed  in  the  desired  direction. 

For  windings  configured  in  a  star  arrangement,  it  can 
be  seen  that  conduction  is  continuous  in  one  leg  while 
commutation  occurs  in  the  other  two  legs.  For  example,  as 
Q1  and  Q5  are  energized  between  30  and  60  degrees  (mechan- 
ical), current  flows  down  leg  A  and  into  leg  B.  At  the  same 
time,  the  current  through  leg  C  due  to  the  energy  storage 
behavior  of  an  inductor  decays  to  zero  through  diode,  D3 
(see  Figure  U.1).  in  the  next  sequence,  60  to  90  degrees, 
C6  is  energized  and  Q5  switched  off.  Current  has  reached  a 
steady  state  flow  through  leg  A,  but  now  it  flows  down  leg  C 
and  the  current  in  leg  B  decays  to  zero  through  D2.  In  the 
next  sequence,  60  to  120  degrees,  Q2  is  switched  on  and  Q 1 
switched  off.  Current  has  reached  steady  state  flow  through 
leg  C  but  now  current  builds  up  and  flows  from  leg  B  while 
the  current  in  leg  A  decays  through  le:j  C.   This  same  action 
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repeats  in  subsequent  commutation  states.  Looking  at  joint 
C  in  Figure  4.1,  it  can  be  seen  that,  since  current  is 
always  continuous  in  one  leg  whether  flowing  into  or  from 
node  D,  Kirchoff s  Current  Law  requires  that  the  net  current 
flowing  in  the  other  two  legs  must  eguai  the  flow  in  the 
continuous  leg.  This  permits  the  following  important 
simplification  which  was  needed  later  in  the  modeling 
process.  Namely,  the  three-legged  stator  can  be  approxi- 
mated by  just  two  windings  serially  configured.  One  leg 
will  always  accurately  reflect  steady  state  flow  conditions 
while  the  other  will  be  treated  as  if  it  were  in  steady 
state  because  the  build  up  in  one  leg  is  balanced  by  decay 
in  the  other.  It  is  recognized  that  this  model  is  a  first 
approximation  and  that  future  modeles  will  require  addi- 
tional refinements  in  order  to  adequately  represent  the 
effect  of  switching  transients  on  the  transistor  and  diode 
elements. 

E.   HALL  EFFECT  SENSCB  FEEDBACK 

As  was  mentioned  in  the  introduction,  the  switching 
action  of  the  commutation  circuitry  is  based  upcn  rotor 
position  feedback  from  Hall  effect  sensors.  In  crder  to 
represent  electronic  switching  and  postional  feedback,  as 
well  as  the  three  phase  winding  configuration,  it  was  neces- 
sary to  modify  the  block  diagram  of  Figure  2.1  to  that  in 
Figure  4.3 

Comparing  Figures  2.1  and  4.3,  it  is  seen  that  the 
chief  difference  is  in  the  switching  logic  block  and  the 
addition  of  two  more  transfer  function  blocks  to  account  for 
the  additional  windirgs.  To  better  understand  the  logic 
relationship  between  the  sensors  and  the  current  switches, 
the  truth  table  for  both  counterclockwise  and  clockwise 
rotation  is  presented  in  Table  II. 
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Figure  4.3    Equivalent  Circuit  With  Electronic  Commutation. 


A  brief  review  of  the  action  af  a  Hall  effect  senscr  is 
in  order  at  this  point.  Consider  a  conducting  material  with 
a  constant  current  flow.  when  there  is  no  external  magnetic 
field  present,  the  current  flows  undisturbed  in  a  straight 
line.  But  when  a  magnetic  field  is  applied,  a  deflecting 
force,   the  Lorentz   force  given  by  II   x  B  (where  I   is  the 
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TABLE  II 
Sensor  and  Snitching  Logic 


RPS  A 


RPS  B 


1 

1 

1 

0 

1 

0 

0 

0 

0 

1 

0 

1 

RPS  C 

0 
0 

1 
1 
1 

0 


?H  A 

high 

high 

open 

low 

low 

open 


COUNTERCLOCKWISE 


0 

0 

1 

0 

1 

0 

1 

1 

0 

1 

0 

1 

CLOCKWISE 


low 

open 

high 

high 

open 

low 


PH  B 

low 

open 

high 

high 

open 

low 


high 

high 

open 

low 

low 

open 


FK  C 

open 

low 

lew 

open 
high 
high 


open 

lew 

low 

open 
high 
h  i  c  h 


current,  1  is  the  length  of  the  conductor  and  B  is  the 
magnetic  field)  causes  the  current  to  bend  and  electrons  to 
pile  up  on  one  side  and  positive  charges,  or  holes,  to  do 
the  same  on  the  other  side.  Thus  a  transverse  Hall  poten- 
tial difference,  Vh,  develops  -  across  the  conductor  (see 
Figure  4.4) . 

Figure  4.5  shows  that  in  the  presence  of  a  sinusoidal 
magnetic  field  produced  by  a  rotating  magnet,  such  as  cne  of 
the  pole-pairs  of  a  multiple-pole  motor,  the  induced  Hall 
voltage  takes  on  an  ideal  square  wave  shape. 

The  logic  circuitry  of  the  sensor  is  configured  in  such 
a  way  as  to  output  a  logic  level  1  for  positive  Hall  volt- 
ages and  a  0  for  negative  voltages.  Briefly  reviewing  the 
action  of  the  motor:  as  the  shaft  turns,  the  rotor  position 
sensors  send  this  information  to  the  switching  logic  block 
which  then  decides  which  phases  are  to  be  made  high,  lew  or 
open  in  accordance  with  Table  II.  Eguivalently,  the 
voltage.  En,  is  then  applied  to  the  two  appropriate  windings 
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Figure  4.4   The  Hall  Effect 

and  zero  voltage  is  applied  to  the  third  winding  in  agree- 
ment with  the  assumption  that  current  is  steady  in  one  leg 
while  the  commutation  in  the  other  two  legs  is  equivalent  to 
a  continuous  current  in  one  leg.  As  before,  each  voltage 
then  produces  a   current.    Because  this  is  assumed   to  be  a 
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Figure  4.5    Sensor  Logic  Based  on  a  Hall  Effect  Sensor 

linear  system,  the  principle  of  superposition  is  applicable 
and,  thus,  the  torques  developed  by  current  flow  through 
separate  windings  car  be  considered  independently  and  then 
summed  to  give  a  total  motor  torgue.  Beyond  this  point,  the 
brushless  motor  components  are  identical  to  its  brush-type 
counterpart  and  behave  similarly. 
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C.   MCDEI  BEVISION 

CSMP  modeling  of  this  version  of  the  brushless  dc  motor 
called  for  two  modifications  to  the  brush-type  model. 
Essentially,  a  CSN5  procedure  function  was  constructed  to 
act  as  a  subroutine  that  embodies  the  logic  of  the  switching 
action.  The  logic  recognizes  that  if  the  motor  shaft  is 
within  seme  degree  range,  for  example  30  to  60  degrees,  then 
it  sets  the  appropriate  rotor  position  sensors  and  energizes 
the  switches  for  the  appropriate  two  windings  and  deener- 
gizes  the  switch  for  the  third.  In  effect,  only  the  forcing 
voltages  are  passed  to  the  main  program,  whereas  switch  and 
sensor  positions  are  available  for  output  (see  Appendix  A) . 
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V.  VARIABLE  FLOX 

A.   AIR  GAP  FLUX 

The  steady  stat€  behavior  of  both  a  brash-type  and  a 
brushless  dc  motor  has  been  studied  under  the  fundamental 
assumption  that  the  tack  emf  generated  is  a  product  of  tne 
shaft  velocity  times  a  constant,  Kb  (see  Equation  2. 2) . 
likewise,  the  motor  torque  is  equated  to  the  product  of  the 
current  times  a  constant,  Kt  (see  Equation  2.12).  But  as 
Equation  2.11  reveals  in  the  case  of  the  torqae  constant, 
Kt,  the  torque  is  in  fact  the  product  of  some  other  constant 
cf  proportionality,  Kt1,  times  the  air  gap  flux,  <|).  Ihe 
same  is  true  for  the  tack  emf  constant,  Kb;  that  is,  Kb  is 
the  product  of  some  proportionality  constant,  Kb1,  times  the 
air  gap  flux,  <{).  In  both  cases,  the  flux,  $,  has  been 
assumed  to  be  constant  which  is  in  agreement  with  the 
general  practice  when . modeling  permanent  magnet  dc  motors. 
In  point  of  fact,  the  flux  is  not  constant  but  rather  is 
distributed  sinusoidally,  generally  being  the  sum  of  a  sinu- 
soid and  one  or  more  of  its  harmonics.  How  the  motor  is 
physically  constructed  determines  this  flux  distribution. 
Two  of  the  chief  factors  are  the  structure  of  the  magnetic 
pole  faces  and,  most  importantly,  how  the  windings  are 
distributed  in  the  armature.  The  latter  includes  such 
factors  as  the  number  of  slots,  the  number  of  coils  per  slot 
and  the  spacing  of  the  slots  [Eef.  5].  The  careful  designer 
distributes  the  windings  so  as  to  minimize  the  effects  of 
harmonics.  Regardless  of  their  constituents,  these 
composite  sinusoids  have  an  average  value,  of  course,  and  so 
in  practice  it  is  this  average  value  that  becomes  the 
constant  flux  term,  (f: . 
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B.   FIOX  AS  AN  AYERAGE  YALOE 

A  logical  next  step  in  the  design  of  a  more  realistic 
model  was  the  determination  of  an  average  value  for  the 
varible  flux,  $av9/  and  the  corresponding  proportionality 
constant,  KK1  such  that  the  ratio  of  (j)avg  and  KK 1  equals  cue 
and  the  relationships  of  equations  2.2  and  2.5  are 
preserved.  The  equations  for  back  emf,  eb  (t)  ,  and  developed 
torque,  tra(t),  now  become 

eb(t)  =  Kbp  *  wm(t)  (e^n  5.1) 

tm  (t)  =  Ktp  *  i(t)  (eqn  5.2) 

where  the  back  emf  constant,  Kbp,  and  the  torque  constant, 
Ktp,  are  now  expressed  as 

Kbp  =  Kb  *  ($avg/KK1)  (eqn  5.3) 

Ktp  =  Kt  *  ((j)avg/KK1)  (eqn  5.4) 

Assuming  as  a  first  approximation  that  the  air  gap  flux 
varies  as  a  simple  sinusoid  of  unit  magnitude,  the  first 
step  toward  finding  the  total  system  flux,  $avg,  was  to 
determine  the  relationship  of  the  individual  flux  patterns 
in  each  of  the  three  phase  windings.  Analysis  of  a  four- 
pole  magnet  rotating  under  three  windings  separated  by  120 
degrees    (electrical)    revealed    the    following    relationships 

(J)a    =    sin(28    +    T7/6)  (eqn    5.5) 

(J)b    =    sin(29    +    9TT/6)  (eqn    5.6) 

<{)c    =    sin  (20   +    5TT/6)  (eqn    5.7) 
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where  (J) a,  $o  and  $c  are  the  flux  contributions  of  the 
respective  phases.  Eeturning  to  an  earlier  assumption  that 
the  circuit  can  be  approximated  as  two  windings  in  series, 
ignoring  the  third,  then  the  average  flux  can  alsc  be 
approximated  from  the  sum  of  the  flux  developed  in  two  wind- 
ings at  a  time.  This  approximation  is  further  supported  by 
the  fact  that  the  flux  in  the  decay  leg  has,  in  fact,  an 
average  value  of  zero  during  every  commutation  interval. 
From  th€  switching  sequence  of  Figure  4.1,  the  two 
conducting  legs  were  identified  for  each  interval  and,  from 
Equations  5.5  through  5.7,  their  respective  flux  values 
calculated  and  summed.  For  every  interval,  the  resulting 
flux  waveshape  was  identical  and  had  an  average  value  of 
(|)avg  =  1.631  (see  Figure  5.1). 

Since  KK1  merely  normalizes  (J)avg,  its  value  was  easily 
determined  to  be  KK1  =  0.613.  The  appropriate  substitutions 
were  then  made  in  the  CSMP  model  (see  Appendix  A  for 
Revision  Two) .  Since  this  amounted  to  nothing  more  than  a 
rearrangement  of  constants,  the  model's  output  behavior 
remained  unchanged.  At  this  stage,  the  model  simulates  the 
behavior  of  a  brushless  dc  motor  being  electronically  ccmmu- 
tated  and  exhibiting  constant,  steady  state  developed  torque 
and  motor  speed  due  to  a  fixed,  average  value  air  gap  flux. 

C.   SINUSOIDAL  FLUX 

The  final  stage  of  this  thesis  was  to  examine  the 
steady  state  behavior  of  the  brushless  dc  motor  under  vari- 
able air  gap  flux  conditions.  As  a  first  approximation,  the 
flux  was  modeled  as  a  simple  sinusoid  and  the  flux  relation- 
ships are  those  given  in  equations  5.5  through  5.7  As 
before,  only  the  flux  relating  to  the  two  conducting  legs 
was  considered.  This  composite  variable  flux,  <j)v1,  is  shown 
in  Figure  5.1  Equations  5.3  and  5.4  had  to  be  modified  to 
account  for  the  different  value  of  flux. 
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Figure  5.1    Composite  Flux  Variation  for  Two  Windings 


Kbpp  =  Kb  *  (<j)v1/KK2) 


Ktpp    =    Kt    *     (<{>v1/KK2) 


(egn  5.8) 
(egn  5.9) 


Recall  that  the  objective  was  to  keep  the  values  of  Kbpp  and 
Ktpp  as  near  as  possible  to  those  of  Kb  and  Kt,  respec- 
tively. In  order  to  roughly  normalize  (f) v  1  #  the  constant 
value  was  determined  as  KK2  =  1.655.  Since  (j)v  1  was  now  a 
variable  guantity,  it  should  be  recognized  that  Kbpp  and 
Ktpp  were  no  longer  strictly  constants  and  were  better 
expressed  as  the  following  functions. 


Kbpp  =  Kb  ($v1) 


(egn  5.  10) 
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Ktpp  =  Kt  (<pv  1)  (egn  5.  11) 

The  equations  for   back  emf ,   eb  (t) ,      and   developed  torque, 
tm { t)  ,  now  became 

eb(t,(J)v1)  =  Kbpp  *  wra(t)  (egn  5.12) 

tm(t,(J)v1)  =  Ktpp  *  i(t)  (eqn  5.13) 

Eecause  the  total  flux  was  no  longer  a  constant, 
average  value,  it  was  expected  that  variation  in  developed 
torque  (and,  conseguently,  in  motor  speed)  would  result. 
Furthermore,  as  the  system  block  diagram  of  Figure  2.1  indi- 
cates, torgue  is  converted  by  the  motor  transfer  function 
into  a  rotational  speed.  This  is  a  process  that  is  equiva- 
lent to  inertia  filtering  of  the  torque  ripple  that  results 
from  variable  flux.  For  the  purposes  of  this  thesis,  ripjle 
is  defined  as  the  ratio  of  the  amount  of  variation  from  the 
maximum  to  the  maximum  itself.  Having  made  the  appropriate 
substitutions  in  the  model  (see  Appendix  A  for  Revision 
Three)  ,  the  program  was  run  under  .no-load  conditions.  It 
was  observed  that  there  was,  in  fact,  ripple  in  both  the 
developed  torgue  and  motor  speed  and  that  some  inertia 
filtering  did  occur.  The  next  logical  step  was  to  lead  the 
motor  and  study  the  effect  on  torgue  and  speed  ripple. 
Table  III  indicates  that  under  no-load  conditions  where 
speed  was  highest,  the  percentage  of  ripple  in  the  output 
speed  was  at  a  minimum  but  that  torque  ripple  was  at  a 
maximum. 

It  can  be  seen  that  as  the  load  increased  the  amount  of 
ripple  in  the  developed  torque  decreased  from  about  8  4?  to 
8%.  A  closer  examination  of  the  torque  behavior  under 
varying  loads  revealed  that,  even  as  the  motor  torque  neces- 
sarily grew  with  increasing  loads,    in  every  case  there  was 
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TABLE  III 
Torque  and  Speed  Ripple  Due  to  Sinusoidal  Flux 


MOTOR  MOTOR 

10AD  (oz-in)  SPEED 

0.0  2557 

32.0  2087 

€4.0  1617 

96.0  1147 


%    RIPPLE 
(EFM)       TORQUE"    "SPEED 


33.9 

1.  1 

21.  1 

1.6 

12.2 

2.6 

7.8 

4.6 

approximately  a  constant  8-9  oz-in  variation.  This  varia- 
tion can  be  attributed  primarily  to  the  flux  variation  in 
Ktpp  and  to  a  lesser  degree  to  the  variation  in  i  (t)  caused 
ty  the  variation  in  tack  emf.  Thus  the  magnitude  of  torque 
variation  remained  about  the  same  for  all  loads  though  the 
percentage  ripple  appears  to  indicate  otherwise.  In  the 
case  cf  motor  speed,  the  amount  of  ripple  appeared  to 
increase  slightly  with  increasing  load  and  thus  decreasing 
speed.  In  this  instance,  the  magnitude  of  variation  did 
increase  a  small  amount  with  decreasing  speed  -  from  28  rpm 
at  no-load  to  52  rpm  at  96  oz-ins  -  and  suggests  that  speed 
ripple,  at  least,  is  somewhat  dependent  on  load.  This  is 
due  in  part  to  the  fact  that  as  the  motor  slows,  there  is 
less  rotor  momentum  and  thus  less  inertia  filtering.  In 
addition,  at  slower  speeds  there  are  less  frictional  effects 
which  at  higher  speeds  also  tend  to  filter  out  ripple.  In 
general,  it  makes  intuitive  sense  that  at  higher  speeds 
frequency  variations  are  harder  to  discern  than  at  slower 
speeds. 

D.   HARHCNIC  FLUX 

Having  observed  the  effects  of  simple  sinusoidal  flux 
variation,  the  next  step  was  to  replace  the  single  sinusoid 
with  a  flux  composed  of  the  sum  of  a  sinusoid  and  its 
harmonics  which  is   mere  representative  of  the   air  gap  flux 
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patterns  in  an  actual  motor.  The  back  emf  of  a  typical 
commercial  brushless  dc  motor  is  given  in  Figure  5.2  where 
the  voltage  between  two  terminals  was  measured  with  the 
motor  rotating  at  1200  rpm. 

A  harmonics  analysis  of  this  waveshape  shows  that  the 
principle  harmonic  was  the  fifth  and  that  the  waveshape  was 
closely  described  by  the  expression 

3.0  sin  (9)  +  0.59  sin(50)  (eqn  5.  14) 

tfhen  substituted  into  equations  5.5  thru  5.7,  and  the 
different  angular  speed  taken  into  account,  the  following 
relationships    for  flux   result 

(()a    =    3.0*sin  (2©+ir/6)    +    0.  59*sin  (1 09  +  5 1T/6)  (egn    5.15) 

$b    =    3.0*sin  (29+9  7T/6)     +    0.  59*sin  (1  00  +  9  V/6)  (egn    5.16) 

(J)c    =    3.0*sin  {20+5  7T/6)    +    0.  59*sin  ( 1  00+  V/6)  (egn    5.17) 

The  same  procedures  used  in  the  simple  sinusoid  model 
were  then  applied  to  the  harmonics  case  where  the  composite 
flux  changes  to  $v2,  a  flux  factor  that  varied  to  a  slightly 
greater  degree  than  did  <j)v1.  Equations  5.3  and  5.9  had  to 
be   modified   as    follows 

Kbppp    =   Kb    *     ((f)v2/KK3)  (egn    5.  13) 

Ktppp    =    Kt    *     ((})V2/KK3)  (egn    5.  19) 

where  KK3  approximately  normalizes  the  flux  and  had  the 
value  KK3  =  5.38.  As  before,  Kbppp  and  Ktppp  varied  with 
the  flux  and  so  the  new  equations  for  back  emf  and  developed 
torgue    became 
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Figure  5.2    BacX  EMF  laveform  Due  to  Harmonic  Flux 


eb(t,<j)v2)  =  Kb  ($v2)  *  wm(t) 


eb(t,$v2)  =  Kbppp  *  wra  (t) 


(eqn  5.20) 
(egn  5. 21) 
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131(1,(1)72)  =  Kt  (<t)v2)  *  i  (t)  (eqn  5.22) 

tm(t,(})v2)  =  Ktppp  *  i  (t)  (egn  5.23) 

The  appropriate  substitutions  were  once  more  made  in  the 
model  (see  Appendix  A  for  Revision  Four)  and  the  model  was 
run  under  various  loads.  As  anticipated,  the  variable  flux, 
<£v2,  produced  ripple  in  the  developed  "torque  and  motor 
speed.  Table  IV  shews  the  resultant  values  for  various 
loads . 

TABLE  IY 
Torque  and  Speed  Ripple  Due  to  Harmonic  Flux 


MOTOR  MOTOR  %  RIPPLE 

IOAD  (oz-in)  SPEED  (rpm)       TORQUE      SPEED 

0.0  2557  94.6        2. a 

32.0  2087  33.5        3.7 

64.0  1617  22.5        5.6 

96.0  1147  14.7        9.8 


Comparing  Tables  III  and  IV,  it  can  be  seen  that  the  values 
for  torque  and  speed  ripple  are  of  similar  magnitudes.  The 
somewhat  higher  values  for  torque  are  due  to  the  fact  that 
in  the  (J)v2  case  there  was  a  16-21  Oz-in  variation  as  opposed 
to  an  average  8-9  cz-in  variation  in  the  <f)v1  case.  The 
speed  ripple  behavior  was  much  like  the  (J)v1  case  with 
proportionate  increases  resulting  from  the  increased  torque 
ripple.  Thus,  for  both  torque  and  speed  ripple  due  to  (f>v2, 
the  model  performed  as  expected;  that  is,  slightly  higher 
values  of  ripple  occurred  which  were  clearly  attributable  to 
the  increased  variability  of  the  flux  pattern  of  composite 
sinusoidal  harmonics. 

An  important  observation   made  during  runs  of   the  tyv} 
model  came  as  a  conseguence   of  the  motor  speed  variability, 
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particularly    when      the    motor      was    unloaded.  Since    no-load 

speed  was  a  function  of  Kb  as  given  in  equation  3.4,  the 
variatle  Kbpp  had  the  effect  of  increasing  the  upper-end 
speed  beyond  the  2557  rpm  baseline  given  in  Figure  3.1  This 
increased  magnitude  ccupled  with  the  swing  of  speed  values 
created  an  undesirable  effect  in  the  unloaded  conditior.  In 
effect,  at  certain  points  in  time  the  additional  speed 
produced  back  emf  values  that  exceeded  the  input  voltage 
which  in  turn  resulted  in  negative  current  and  developed 
torque  values.  Though  these  were  transient  negative  values 
and  had  small  effect  on  the  steady  state  operation  of  the 
motor,  still  their  presence  was  obviously  unrealistic.  The 
solution  to  returning  the  motor  to  all  positive  values  of 
current  and  torque  was  to  increase  the  friction  filtering 
effect  by  adjustment  of  the  viscous  friction  coefficient, 
Bm.  It  will  be  recalled  from  Chapter  Three  that  a  value  for 
the  Bm  of  the  commercial  motor  being  modelled  was  not  avail- 
able and  that  for  the  basic  model  a  value  for  Bm  was  derived 
from  the  curve  fitting  process.  Thus,  it  seemed  reasonable 
and    permissable    to   repeat    the    procedure    here. 

Since  it  makes  sense  that  increasing  the  friction  within 
a  system  will  slew  the  system  down,  in  this  case  increasing 
the  value  of  Bm  did  just  that.  By  increasing  Em  from 
0.000  15  to  0.045  oz-in/rad/s,  the  upper-end  of  the  motor 
speed  was  reduced  to  a  value  that  eliminated  the  undesirable 
negative  values.  This  caused  the  no-load  speed  to  be 
lowered  slightly  belcw  the  desired  value  of  2557  rpms. 
Recalling  that  no-lcad  speed  is  fixed  by  the  back  emf 
constant,  Kb,  this  parameter  was  adjusted  from  0.1120  to 
0.1089  vclt/rad/s  to  bring  the  back  up  to  speed.  Since  this 
adjustment  was  less  than  3%  and  the  manufacturer  specifica- 
tions allowed  a  105?  +/-  measurement  error  margin,  this 
adjustment    to      Kb    was  permissable.  Of   course,         a    propor- 

tionate   adjustment    to  the    torque    constant,    Kt,    was   also   made 
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since  Kt  and  Kb  must  maintain  a  constant  relationship  as 
explained  in  Chapter  Three.  The  effect  on  motor  speed 
ripple  was  neglible.  Reviewing  the  ({)v2  model  at  this  joint 
revealed  that  similar  behavior  in  the  no-load  state  occurred 
but  en  a  larger  scale  due  to  the  greater  variability  of  <f>v2. 
A  sinilar  set  of  adjustments  was  necessary  to  eliminate  the 
negative  current  and  torgue  values  in  the  unloaded  condi- 
tion. It  reguired  an  adjustment  of  Bm  from  0.00015  to  0.045 
oz-in/rad/s  in  order  to  ensure  that  all  values  were  posi- 
tive. In  addition,  the  back  emf  constant,  Kb,  was  adjusted 
to  0.1192  volt/rad/s  and  a  proportionate  adjustment  to  Kt 
was  made.   Again,  all  adjustments  were  within  tolerances. 

Of  course,  it  was  earlier  established  in  Eguation  3.7 
that  changes  to  Kb  and  Kt  would  change  the  slope  of  the 
speed-torgue  curve.  Thus,  in  both  cases  the  slope  was 
returned  to  its  original  value  by  adjusting  the  motor  resis- 
tance, R.  In  the  case  of  sinusoidal  variable  flux,  E  was 
reduced  from  2.74  to  2.70  ohms.  For  harmonic  variable  flux, 
R  was  increased  from  2.74  to  3.0  ohms.  Again,  both  adjust- 
ments were  within  the  10%  ♦/-  allowable  deviation.  It  was 
also  observed  that  since  both  Kb  and  Kt  no  longer  remained 
constant  but  rather  varied  within  small  ranges,  the  slope 
was  not  constant  and  varied  accordingly  in  both  the  ij)v1  and 
$v2  models  with  the  greatest  change  occurring  when  Kt  and  Kt 
simultaneously  reached  the  minimum  values  of  their  respec- 
tive ranges.  Bearing  in  mind  that  these  deviations  in  slope 
are  fairly  transient,  their  average  values  are  most  impor- 
tant and  are  the  approximate  values  from  which  the  final 
adjustments  to  R  were  based. 
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VI.  SOMHARI 

A.   REMARKS  AND  CONC1DSIONS 

CSMP  language  is  a  convenient  tool  for  modeling  a 
fcrushless  DC  motor.  Positional  sensor  feedback  and 
switching  logic  being  functions  of  time,  the  action  cf  elec- 
tronic commutation  can  be  nicely  simulated  in  CSMP.  The 
model  can  be  fairly  easily  modified  in  order  to  study  ether 
motor  configurations.  For  example,  the  three-phase  star 
configuration  used  in  this  thesis  could  easily  be  changed  to 
a  grounded  neutral,  point  D  in  figure  4.1,  and  the  same 
types  of  analysis  performed.  In  future  studies  where  the 
complexity  of  the  mcdel  is  increased  in  order  to  include 
more  design  detail,  the  advantage  of  CSMP  will  become  even 
more  apparent. 

In  the  course  of  simulating  the  effects  of  different 
air  gap  flux  variations,  it  appeared  necessary  to  make  minor 
parameter  adjustments  in  order  to  offset  the  feedback 
effects  cf  torgue  and  speed  ripple.  Within  the  scope  of 
this  thesis,  adjusting  the  steady  state  performance  was 
reasonable  and  justifiable.  In  a  larger  context,  however, 
various  control  systems,  particularly  torque  generation 
controllers,  are  specifically  designed  to  reduce  tcrgue 
ripple  and  its  effects.  So,  in  at  least  one  sense,  these 
adjustments  were  somewhat  artificial.  To  the  extent  that 
parameter  adjustment  is  necessary  and  desirable,  it  must  be 
understood  that  this  is  not  a  simple  procedure  since  adjust- 
ment of  one  parameter  almost  always  requires  adjustment  of 
one  or  more  others.  These  adjustments  in  turn  necessitate 
readjustment  of  the  original  parameter.  The  process  is  thus 
an  iterative  one.   This  model  permits  this  kind  of  parameter 
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adjustment,   but  it  is  a  long  and  tedious  process  that  would 
be  tetter  performed  by  some  sort  of  optimization  algorithm. 

E.   RECCflMENDATIONS  ECR  FUTURE  STUDY 

The  simulation  cf  ripple  effects  is  particulary  rele- 
vant to  the  control  engineer  who  must  design  the  control 
circuitry  and  power  conditioner  to  meet  specific  performance 
reguirements.  This  points  to  several  areas  for  future 
study.  One  area  is  the  control  of  torgue  generation  ana  the 
several  principle  configurations  used  in  brushless  DC 
motors.  Two  commonly  used  methods  are  the  sinusoidal  torgue 
generation  principle  and  the  trapezoidal  torgue  scheme. 
[Ref.  5].  In  addition  to  torgue  generation  control,  a 
second  area  for  investigation  is  power  control  of  brushless 
EC  motors.  In  this  thesis,  output  power  was  unregulated  and 
was  studied  simply  as  a  function  of  a  constant  input 
voltage.  In  practice,  however,  power  is  produced  and 
controlled  by  varying  the  supply  voltage.  Various  schemes 
exist  for  power  control;  amoung  these  are  linear  transistor 
control  and  pulse-width  or  pulse-f reguency  control. 
[Ref.  5]. 

A  very  important  area  that  this  thesis  did  not  address 
was  that  of  the  switching  transients  that  occur  as  a  conse- 
guence  of  commutation.  During  each  switching  interval,  each 
winding  is  an  inductor  that  stores  energy  which  causes  a 
significant  amount  of  current  to  flow  at  the  instant  of 
switching.  This  can  have  serious  effects,  especially  if 
breakdown  conditions  exist  amoung  the  commutation  transis- 
tors and  other  control  circuit  elements.  In  this  same  area, 
switching  transients  also  occur  within  the  control  circuitry 
itself.  Thus,  a  thorough  investigation  and  understanding  of. 
the  behavior  of  power  transistors,  diodes  and  other  control 
elements  used  with  electronic  commutation  needs  to  be  done. 
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APPENDIX  A 
LISTING  OF  MODEL  PROGRAMS 

A.   BASIC  PROTOTYPE  MODEL 

The  program  in  this  section  is  a  basic  program  that 
simulates  the  action  of  a  standard  brush-type  DC  motor.  It 
is  the  prototype  from  which  the  brushless  DC  motor  model  and 
its  several  revisions  are  derived. 
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B.   REVISION  ONE 

This  is  the  first  revision  of  the  basic  prototype.  This 
model  simulates  the  switching  logic  of  a  brushiess  EC  motor 
which  is  based  on  feedback  from  Hall  effect  sensors.  In 
addition,  the  windings  are  treated  independently  and  their 
contributions  to  developed  torgue  are  superposed.  Since  the 
superposition  results  in  twice  as  much  current  flew  as  a 
single  lumped  coil,  the  total  current  is  halved  to  retain 
the  same  overall  motor  behavior  as  in  the  prototype. 
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C.   REVISION  TWO 

In  the  following  program,  the  motor's  air  gap  flux  is 
treated  as  the  average  value  of  the  sum  of  the  flux  acting 
in  two  windings  as  explained  in  Chapter  Five.  For  conven- 
ience of  study,  both  windings  are  again  lumped  into  a  single 
winding  as  in  the  prototype.  This  facilitates  a  closer 
focussing  on  the  effects  of  variable  flux  on  motor  tehavicr. 
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D.   REVISION  THREE 

The  purpose  of  this  version  of   the  program  is  treat  the 

air  cap  flux  as  varying   in  simple  sinusoidal  fashion.  The 

total  system   flux  is   treated  as  the   algebraic  sum   of  the 
flux  developed  in  the  two  active  windings. 
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E.   EEVISICH  FOUR 

This  version  of  the  model  treats  the  flux  as  varying 
according  to  the  sura  cf  a  sinusoid  of  fundamental  frequency 
and  its  fifth  harmonic  as  explained  in  Chapter  Five.  The 
total  flux  is  again  approximated  as  the  algebraic  sue  cf  the 
flux  developed  in  two  windings  at  a  time. 
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A£EMJ2II  b 

SAMPLE  OOTPOT 

The  following  is  a  sample  of  the  type  of  information 
the  mcdel  can  provide.  In  this  particular  case,  the  fourth 
revision  of  the  basic  prototype,  the  harmonic  flux  case,  was 
used.  For  completeness,  the  motor  was  run  at  loads  fron;  0.0 
up  through  96.0  oz-in,  which  was  the  range  of  leads  with 
which  all  analysis  was  performed.  In  addition  to  such  motor 
variables  as  speed,  current  and  power,  a  sampling  of  the 
feedback  from  Hall  effect  sensors  as  well  as  the  switching 
action  of  the  power  transistors  is  included. 
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